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Abstract 

We study the bulk viscosity of strange quark matter (SQM) in a quasiparti- 
cle model at finite chemical potential by extrapolating the previous quasiparticle 
model of finite temperature lattice QCD. The more proper bulk viscosity coeffi- 
cient can be given in this model where chemical potential fi and coupling constant 
g are interdependent. We also apply our result to determine the critical rotation of 
strange stars by r-mode instability window. Our model is compatible to the mil- 
lisecond pulsar data for a wide range of mass and radius of the stars. 
PACS numbers: 97.60.Jd, 12.38.Mh, 97.60.Gb 

1 INTRODUCTION 

According to Witten's conjecture that the strange quark matter (SQM) composed 
of roughly equal numbers of up, down and strange quarks might be absolutely stable or 
metastable phase of nuclear matter in 1984 1 1 1, strange stars might exist in the universe, 
and their structure and properties have been widely studied j2l|3]|4l|5]|6|El- 

How to distinguish strange stars from neutron stars has been one of the important 
issues. Madsen pointed out that the r-mode instability in the relativistic stars at all rotat- 
ing rates may provide a distinguishable signature 1 8 1. The main reason is that the bulk 
viscosity coefficient of SQM is larger than that of neutron matter, and hence suppresses 
the r-modes. Since Wang and Lu |9| found that the non-leptonic weak reaction domi- 
nates the bulk viscosity of SQM, some investigations have tried to calculate the relevant 
viscosity coefficient of SQM. The equation of state (EOS) of strange stars has also been 
studied as the base of studying the viscosity. Sawyer 1 2 1 and Madsen 1 3 1 completed the 
calculations of bulk viscosity for the ideal quark gas. We had investigated the viscosity 
of interacting SQM in quasiparticle description by regarding the coupling constant as 
an independent parameter and found the improved result is of importance for astro- 
physical relevance 1 1 1 . However it is desirable to have a more consistent investigation 
by considering the realistic running coupling constant for SQM in the interior of com- 
pact stars. According to the regulations, the realistic coupling constant and EOS need 
nonperturbative evaluations, i. e. lattice QCD calculations. However, the present lattice 
QCD calculations are yet restricted to zero chemical potential /i. Attempts to extend the 
lattice calculations systematically to non-zero chemical potential are underway. Being 
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aware of the urgent need for the EOS, many works have suggested an approach based 
on a quasiparticle description of quarks and gluons, to map available lattice data from 
li = to finite values of fi and to small temperatures II II . Consequently, the run- 
ning coupling in the asymptotic limit of large chemical potential can be also simulated 
from that of large temperatures at zero chemical potential, which will also help us to 
overcome the difficulty of calculating EOS in non-perturbation regime. We thus can 
consistently evaluate the bulk viscosity in the light of the simulated EOS and running 
coupling constant. We will see below the coupling among quarks influences remark- 
ably the bulk viscosity in SQM. Our theoretical output is appropriate for a wide range 
of stellar parameters. 

We organize this paper as follows. In Sec. 2, we introduce the equation of state and 
formulate the quasiparticle model at finite chemical. In Sec. 3, the bulk viscosity of 
SQM in the quasiparticle model is derived, which arises from the nonleptonic weak 
interaction. In Sec. 4, we probe the application of our model. Finally, we summarize 
our conclusion and discussion in Sec. 5. 



2 THE QUASI PARTICLE MODEL 



Matter in local thermodynamical equilibrium can be described by its EOS which 
represents an important interrelation of state variables. For the EOS of SQM, non- 
perturbative methods as lattice QCD should be applied. As a matter of fact, these 
simulations are presently restricted to vanishing chemical potential /i, and the imple- 
mentation of physical quark masses is still too expensive numerically. So we often 
consider phenomenological model. The bag model EOS is popular, but it is in conflict 
with thermodynamic QCD lattice data. We will apply the quasipaticle model and here 
repeat a simple extrapolation of finite temperature lattice to nonzero chemical potential 
relying on thermodynamic consistency II II . 

Asymptotically, the collective behavior of the plasma can be interpreted in terms of 
quasiparticle excitations with a dispersion relation ujf (k) « nij+k 2 andmf = m^+II* 
depending on the rest mass and the leading order on-shell self-energies 11 21 . 
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where /i 9 denotes the quark chemical potential, andA^ c =3. The pressure of quasiparticle 
system can be decomposed into the contributions of the quasiparticles and their mean- 
field interaction B, 



p(T,^;m^) = ^pJr^O);™, 2 ] -5(11*), 



(3) 



where p i — ±diT J d 3 k / (2tt) 3 \n(l ± exp{ — [uji — fii)/T}) are the contributions of 
the gluons (with vanishing chemical potential) and the quarks (for the anti-quarks, the 
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chemical potential differs in the sign). And d g = 2(N 2 — 1) and d q — 2N C count the 
degrees of freedom 1131 . The function B(H*) is determined from a thermodynamical 
self-consistency condition, via 



dB _ dp J (T,^ J ;m 2 j 

an* ~ dm) 



(4) 



Furthermore, the stationarity implies that the entropy and the particle densities are 
given by the sum of the quasiparticle contributions, 
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while the energy density has the form e = £\ <?; + B. By comparison with lattice data, 
at n = 0, Peshier lll II has tested the quasiparticle approach which is an appropriate 
description even close to the confinement transition, with the effective coupling in Eq 
(1) and (2) nonperturbatively parameterized by 

S 2 (7> = 0) = (6) 



(lW c -2tf/)]n(^y 



interpolating to the asymptotic limit of QCD. 

Encouraged by the successful quasiparticle description of the \i = lattice data, 
the model is now extrapolated to finite chemical potential. In general, the pressure is a 
function of the state variables T and /i. As a direct consequence thereof, the Maxwell 
relation implies for the quasiparticle model 
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which is the integrability condition for the function B defined by Eq (4). Following 
directly from principles of thermodynamics, we can get a flow equation for the effec- 
tive coupling with II* depending on g 2 . This flow equation is a quasilinear partial 
differential equation of the form 

dg 2 dg 2 

with the coefficients and b depending on T, /i and g 2 . 

The flow of the effective coupling is elliptic in the nonperturbative regime, thus 
mapping the /i = axis, where g 2 can be determined from lattice data, into the /it — T 
plane. As an example, the nonperturbative flow of the coupling of the Nf = 2 + 1 
system is shown in Figure. 1. Here we use the parameters A = 6.6 and T s = —0.78T C 
1141 . Furthermore, the running coupling constant as a function of chemical potential 
g(T = 0, fi) is displayed in Figure. 2. For convenience later, we can approximately 
formulate g(T = 0, /i) as 

g 2 (T=0, f ,) = ^ ™— . (9) 

(UN c -ZNf)M& 2 
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Based on the extended coupling constant, the EOS of SQM in the interior of compact 
stars is immediately obtained from Eqs (1), (2), (3) and (4). 
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Figure 1: The characteristics of the coupling flow equation (8) for the QCD plasma 
with Nf = 2 + 1 flavors. 



3 BULK VISCOSITY 

Viscosity is important for describing the transport property of matter and can be 
generally calculated from either quantum field theory or kinetic theory 1 15 1. The bulk 
viscosity of SQM mainly arises from the nonleptonic weak interactions 1 9 1 

u + d< > s + u. (10) 

The importance of dissipation due to the reaction (10) was first stressed by Wang and 
Lu in the case of neutron stars with quark cores. Afterward, a series of investigations 
have tried to calculate the viscosity coefficient of SQM in MIT bag model and showed a 
huge bulk viscosity of strange matter relative to nuclear matter. Furthermore, Zheng et 
al 1 10 1 considered the medium effect on the bulk viscosity in quasiparticle description 
and found the viscosity is few^tens times larger than that of non-interacting quark gas 
due to the small influence of medium on the equation of state, where coupling constant 
g is regarded as an independent parameter of chemical potential. In fact, // and g should 
be interdependent according to the results of lattice calculation. We here reevaluate the 
bulk viscosity coefficient adopting the description in Sec. 2. 
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Figure 2: The coupling constant g(T = 0,/j.) as a function of /i . 



In quasiparticle approximation, we can give the formulae of the bulk viscosity in 
accordance with Eqs (1), (2), (3) and the reactional rate of reaction (10) at finite chem- 
ical potential when temperature is small. Zheng et al 1 10 1 have made derivations of 
bulk viscosity in quasiparticle description. We can apply the result, 
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where k Fi = tf-mtfl*. Q = ^ - ^ - - the 
expression of C, has been reevaluated, differing from that of Zheng et al 1 10 1 because 
g is the function of ijl. The effective mass nii was given in Sec. 2, the current masses 
vanish for up and down quarks. We continue to use the formula of the reactional rate 
adopted by Madsen 1 3 1 
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When the temperature is high enough, i. e. 2nT ^ <5/i, the cubic term <5/i 3 in Eq 
(12) can be neglected. We can obtain an analytical result which is similar to Madsen's 
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expression |3|, 

a*T 2 

C = uj 2 + f3*T 4 (U) 

and 

a* = 9.39 x 10 22 ^ f^N ~ ^) (^cm- 1 ^ 1 ), (15) 

2 

(s- 2 ), (16) 

here w = — , a* and /?* remarkably differ from Madsen's. 

In general, Eqs (1 1), (12) and (13) must be solved numerically due to the existence 
of 6/j, 3 in rate (12). The results of such calculations are shown in Figure 3 and Figure 
4. The magnitude of viscosity coefficient increases with chemical potential for the 
low-temperature case such as T = 1CP 4 MeV in Figure 3. It is just contrary to the 
high-temperature case as T = 1CP 1 MeV in Figure 3. This implies a shift of the 
the maximal viscosity to the low-temperature when chemical potential increases as 
shown in Figure 4. By comparison, the viscosity has more complicated dependence on 
coupling constant than that of ref 1 10) . 
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Figure 3: Bulk viscosity coefficient of different /x and T, mo s =200 MeV, r=0.001 s. 
The lower curves denote T = 10~ 4 MeV, the upper curves denote T = 10 _1 MeV. 
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Figure 4: Bulk viscosity as function of temperature for different jx, mo s =200 MeV, 
r=0.001 s, = !0 -4 - 

4 CRITICAL ROTATION OF STRANGE STARS 

Given the EOS of SQM, we can get radius-mass characteristic relationship of 
strange stars through the Tolman-Oppenheimer-Volkov(TOV) equation 1 16 1 for energy 
densities up to several times the nuclear saturation density and at temperatures less than 
some 10 MeV, the relation e{p) of /3-stable quark matter can be parametrized ifTTI by 
e = 4B+ap, as estimated numerically by our quasi-particle model, here energy density 
e = e u + ed + e s + e e + B and pressure p = p u + p d + p s +p e — B n . For the chosen 
parameters in Sec. 2 and the bag constant B^ 4 — 145 Mev, we can get a = 3.5 and 

B 1 ^ 4 = 152 MeV in our EOS. Furthermore, we can find corresponding values of the 
parameter^ by changing the bag constant Bq. The values of B have a strong impact on 
the star's mass and radius, obtained by integrating the TOV equation. In Figure. 5, the 
mass as a function of the radius of strange stars is displayed for several values of the 
parameter B. The maximal mass and radius are essentially consistent with canonical 

pulsar data for bJ /4 - 145 Mev and B V4 - 152 MeV. 

Now we focus on the r-mode unstable window for this case. Andersson in 1998 
1171 recognized the existence of unstable r-modes in perfect fluid stars at all rates of 
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rotation due to gravitational wave emission. The r-mode unstable regime of the real- 
istic stars, neutron stars as well as strange stars, depends on the competition between 
the gravitational radiation and various dissipation mechanisms. To plot the instabil- 
ity window of r-mode or obtain the critical rotation frequency for given stellar model, 
we need to work out the characteristic timescales,damping and growing timescales of 
r-mode instability. The time scale for gravity wave emission is 1181 

t g =T G (nGp/n 2 ) 3 , (17) 

where, tq is -3.26 s for n = 1 polytropic EOS, p is the mean density of the star. In 
strange stars, the time scales for shear and bulk viscous dissipations can be respectively 
given 

T5=r s (a s /0.1) 5 / 3 T 9 5/3 , (18) 
T B =T BL {TTGp/n 2 )T 9 2 (19) 

for low-T limit, 

r B =T BH (7rGp/n 2 ) 2 T 2 (20) 
for high-riimit, where t s is 5.37 x 10 8 s for n=l, t bl = 2.83 x 10 3 a* _1 pmf 00 
a s = ^— is QCD fine structure constant, T 9 denotes temperature in units of 10 9 K, 
mioo denotes the current mass of strange quark in units of 100 MeV. Obviously, tb 
is determined with the chemical potential p,d- For given mean density 7j, We can find 
appropriate fid, and g via conservation of baryon number, 

n B = + n d + n s ) (21) 

and 

P = (§) n B , (22) 

here, n b reprensents the baryon number density, and «j = -^k 2 Fi is the flavor number 
density, -| = e ^+ e d+e^ s +e e +B a j s tng energ y per baryon. 

We can evaluate the the critical spin frequency as a function of temperature from 
the equation 

— + — + — = 0. (23) 

tg t s t b 

Figure. 6 shows the regions of r-mode (in)stability in spin frequency- temperature^— 
T) plane for a strange star with mass M = 1.4 Mq and radius R = 10 km. In compari- 
sion, the instability window of usual strange stars is also depicted in this figure, which 
is accomplished by Madsen. We find the medium effect narrows the r-mode instability 
window. The lowest limiting frequency of upper contour is 541 Hz (the corresponding 
period is 1.85 msec), which is more close to the most rapidly spinning pulsars with 
frequencies of 642 Hz and 622 Hz (the periods are 1.56 and 1.61msec). This implies 
a strange star would slow down by gravitational window and spin around in 1.85 msec 
instead of the 2.5 ~ 3 msec expected by Madsen. Figure. 7 describes the instability 
windows for a wide parameters a, B. The dotted curves denote a strange star with 
mass M = 1.0 Mq and radius R = 6.6 km, and the solid curves express a strange star 
with mass M = 2.06 M Q and radius R = 13.5 km. 
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5 SUMMARY AND DISCUSSION 



In accordance with the quasiparticle description, we rederive and reevaluate the 
bulk viscous coefficient of SQM basing on referring // dependence of the coupling 
constant. We find that below few times 10 9 K, the medium effect is enhancing the 
viscosity. The dissipation of fluctuations in dense matter is more effective for this 
situation. We apply our result to study the r-mode instability window of strange stars. 
We find the r-mode instability window is narrowed and the lowest limiting frequency 
is closer to the millisecond pulsars relative to Madsen's one for stars with mass M = 
1.4 M & and radius R — 10 km. Our model also allow for wide frequency changes from 
400 Hz to 900 Hz. The scenario is compatible with pulsar data for a wide range of 
mass (1.0 ~ 2.8 M Q ) and radius (6 ~ 15 km). 
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Figure 5: The dependence of the mass of strange stars on the radius for a = 3.5 and 
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Figure 6: Critical spin frequencies for strange stars as functions of temperature with 
M = I.AMq and R = 10km. The dotted contour stand for Madsen's model, the solid 
curves display our result. 
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Figure 7: Critical spin frequencies for strange stars as functions of temperature. The 
dotted contour stand for a strange star with M = 1.0M Q and R = 6.6km, the solid 
curves display a strange star with M = 2.06M Q and R = 13.5km. 
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